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Abstract. A model of the passive vector field advected by the uncorrelated in time 
Gaussian velocity with power-like covariance is studied by means of the renormalization 
group and the operator product expansion. The structure functions of the admixture 
demonstrate essential power-like dependence on the external scale in the inertial range 
(the case of an anomalous scaling). The method of finding of independent tensor 
invariants in the cases of two and three dimensions is proposed to eliminate linear 
dependencies between the operators entering into the operator product expansions of 
the structure functions. The constructed operator bases, which include the powers of 
the dissipation operator and the enstrophy operator, provide the possibility to calculate 
the exponents of the anomalous scaling. 



Introduction 

The term "anomalous scaling" refers to deviations from the predictions of the 
Kolmogorov theory. Such deviations take place, in particular, in the behaviour of the 
structure functions S n (r) of the turbulent velocity field v in the inertial range I <^ r <ti L, 
where L and I are the external scale and the dissipation length. The structure functions 
of the vector field (p are defined as 



For the field v the Kolmogorov theory predicts S n (r) ~ r n / 3 , while the experiments 
indicate that S n (r) ~ r ™/ 3 -€« with certain nontrivial "anomalous exponents" £ n . 

The Kolmogorov theory states that the only dimensional parameter which 
determines the statistics of the turbulent velocity pulsations is the average rate of 
energy dissipation per unit mass e = z/($d is ), where v is the kinematic viscosity 
and $dis = (diVj + djVi) 2 /2 is the operator of local energy dissipation. In stationary 
state, it is held that e = W with the power of energy pumping W. In this case the 
Kolmogorov hypothesis can be formulated for the structure functions (JTJ) as follows. 
Dimensional analysis leads to the relation S n (r, u, W, L) = (Wr) n ^ 3 R n (r /l,r / X), where 
/ = (z/ 3 /W) 1 ^. Then the Kolmogorov hypothesis states that function R n has finite limit 
in the inertial interval: R n (oo, 0) = const ^ 0. 





(1) 
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At present the possibility of taking limit in the R n (r/l,r/L) on the first argument 
is universally recognized, i.e. the quantity R n {r / L) = R n (oo,r / L) is recognized finite. 
Hence the structure functions are considered to be independent on the viscosity in the IR 
region r>L If the Kolmogorov hypothesis about finiteness of the R n {r/L) in the limit 
r/L —>■ is violated and the asymptotic behaviour is power-like R n (r/L) ~ (r/L)~^ n 
then the anomalous scaling appears. In phenomenological generalizations of the 
Kolmogorov theory the anomalous scaling is usually considered as a result of fluctuations 
of the energy dissipation rate [H Ej . 

The analogy with the theory of critical phenomena suggests itself when the problem 
is formulated as the field-theoretical model. Then the direct analogy to the theory 
of critical phenomena provides possibility to apply the powerful UV-renormalization 
technique to the problem. In this way the first Kolmogorov hypothesis can be proved 
(independence of the structure functions on the viscosity in the inertial range). The 
behaviour of R n (r/L) in the range r/L <C 1 is not determined by the renormalization 
group itself, the method to find it is the operator product expansion (OPE) which gives 
the asymptotic expansion 

5 n (r)ocr"/ 3 j: F (r/L) A M F , (2) 

where summation runs over all possible scalar operator products F (constructed from 
local products of the fields v i(t, x) and their derivatives), A F are the critical dimensions 
of the operators More precisely, Ap are the eigenvalues of the matrix of critical 
dimensions, and summation in (j2J) runs over the eigenvectors of the matrix. The 
operators entering into the OPE are those which appear in the Taylor expansion and 
all the operators that admix to them in renormalization. 

If for each F holds A^ > (as in the theory of phase transitions) then the terms 
in (J2J) determines corrections to the Kolmogorov scaling. If there is an operator in 
(J2J) with A^ < - - "dangerous operator" then the limit r/L — ► in (J2J does not 
exist which leads to the anomalous scaling. Realization of the described scenario meets 
however two difficulties: technical difficulty of calculation of the critical dimensions 
A p and fundamental one — if in the theory a dangerous operator exists then it can be 
shown that with necessity there are infinitely many dangerous operators, the spectrum of 
their critical dimensions is unbounded below (we can not point at "the most dangerous" 
operator). Thus the expansion (J2J) can be useful when summation of the series is possible, 
or when the (J2J contains only finite number of terms due to the model features. 

1. A passive scalar admixture 

Recently, significant progress in the description of the anomalous scaling has been 
achieved in related problems of the turbulent advection of a passive admixture. The 
experiments and computer simulation data demonstrate that the anomalous scaling 
appears not only in the turbulent pulsations of the velocity, but much more it reveals in 
the properties of the field transferred by the turbulent flow 9(t, x) which may be the field 
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of the admixture concentration or the temperature field. The passive scalar advected 
by the turbulent velocity field Vi(t, x) is described by equation 

d t e + (v 3 d 3 )e = u d 2 e + /, (3) 

where v$ is the diflusivity (or the thermal diffusivity) , / is the source of the field 9. 
Significant progress in the description of the anomalous scaling has been achieved in 
a relatively simple model, due to Kraichnan [HE], of a passive scalar advection which 
assumes the velocity field to be delta-correlated in time and Gaussian with covariance 

r dk 

(v i (t + r,x + T)v j (t,x)) = D S{r) J — ^f^(k)iV(&)expikr, (4) 

where -Pjj(k) = 5ij — kikj/k 2 is the transverse projector (the consequence of transversality 
of the velocity), the function N(k) in (^-dimensional space is modelled with power-like 
expression N(k) = k~ d ~ £ and is supposed to be somehow IR regularized on the scale 
k ~ L^ 1 . Here < e < 2 is a kind of Holder exponent which measures the "roughness" 
of the velocity field. In the RG approach, it plays the same role as the parameter 4 — d 
in the theory of critical phenomena. The physical value of parameter e is e = 4/3, for 
which the amplitude factor Do has the same dimension as the energy pumping power 
per unit mass. 

In the model ©-(III), the existence of the anomalous scaling was substantiated and 
the exponents £ n of the structure functions were calculated [HI 13 El E] • The first results 
were obtained using the Hopf equations on the distribution function of the equal time 
fluctuations of the admixture field. The linearity of the equation © for 9 results in 
closed equations (instead of chain of equations) for the equal time correlation functions. 
The analysis of so-called "zero modes" allows to substantiate the anomalous scaling and 
to calculate the anomalous exponents. The pair function is founded exactly and has no 
anomalous scaling (£2 = 0). The exponents of the high-order functions were calculated 
approximately with small parameters e E] and 1/d E] : 

^ n = n(n~2)e/2(d + 2)+0(s 2 ), (5) 

£ n = n (n-2)e/2d + 0(l/d 2 ). (6) 

The analysis of the zero modes gives also interesting results in another case — in the 
Batchelor limit e = 2 PH [TTJ H2 EH . 

Until now the technical difficulties do not allow to calculate higher-order terms of 
the expansions (0) using analysis of the zero modes. The methods of RG and OPE turns 
out to be more effective to find the asymptotic approximation on small e. The terms 
~ e 2 [3 a and ~ e 3 [H] were calculated using these methods. It is essential that in the 
Kraichnan model the finite number of terms contribute to the OPE (j2j) of the structure 
function S n , precisely, those are the powers of the dissipation operator of admixture 
field $*j is with 1 < k < n/2. Joining of asymptotic expansion of the exponents £ n on 
small e (taking three terms of the expansion) with asymptotic approximation near the 
Batchelor limit (e = 2) results in interpolation for £ n (e) which coincides with the data of 
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numerical experiment (within the experimental errors) in the range < e < 2, including 
the physical value e = 4/3 [T4*] . 

The RG method demonstrates the universality of the approach besides technical 
efficiency Using the RG the anomalous exponents of order e 2 in the model of 
compressible fluid were calculated [THj, the model of advection by the Gaussian velocity 
field with finite correlation time was analyzed in [TS] ITTj . the Kazantsev model of the 
magnetic hydrodynamics — in |15j . 

2. A passive vector admixture 

In the present paper we concentrate on the model of advection of a passive vector 
admixture proposed in [TBI EH 1201 E] • The model is the direct generalization of (J3J) and 
is described by the equation 

dm + {vjdj)ipi = is d 2 ipi - diP + f u (7) 

which come out from equation (jHJ) after replacing of the scalar field 9 with the transversal 
vector field y?j and adding of the pressure gradient P to the left hand side to ensure 
transversality diipi = 0. The random force /, in (JZJ) obeys the Gaussian distribution 
with zero mean and the covariance 

(fi(t + r, x + r)fj(t, x)) = 6(t) C^r/L). (8) 

The exact form of the function Cij(r/L) is not important. It only models the 
energy influx resulting from the interaction with large vortices, which compensates the 
dissipation losses. 

Note that all the odd-order structure functions vanish because the equation ((7|) 
is invariant with respect to the substitution ip — > —ip, f — > —f and the fact that / 
is Gaussian (the pressure term can be eliminated from (|7J) by applying the transverse 
projector). 

The model (jlj), ((Zj), (jEJ) can be considered as a rough approximation describing the 
turbulent velocity field, when (p is regarded as the "hard" component of the velocity 
field and v is regarded as the "soft" one. Then in the Navier-Stokes (NS) equation 
for the hard component only the convective transfer by the soft component is taken 
into account and the soft component obeys Gaussian distribution with covariance 
(J3J). These assumptions do not violate Galilean invariance of the model. Since the 
structure functions are Galilean- invariant, it follows that OPE (J2J) contains only 
contributions from the operators which are also Galilean- invariant, just as in the 
stochastic hydrodynamics which is described by the NS equation. 

The model (Jlj), ((Zj), (jEJ) was analyzed by means of the Hopf equation for the pair 
correlation function [T^]. The function turns out to be non-anomalous (£2 = 0) in 
the isotropic case. The model was considered using the methods of RG and OPE 
in [IB]. It was shown, that the main contribution to the OPE of the S^n gives a 
family of scalar operators of the form (<9y?) 2n . Renormalization of the operators is not 
multiplicative, but includes operator mixing. The mixing complicates greatly calculation 
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of the anomalous exponents of the operators; obtaining of an analytic expression of the 
anomalous exponents for arbitrary n do not succeed even in one-loop approximation. In 
|18j . the one-loop approximation was considered for the family of the S4; the exponents 
were evaluated in the linear on e approximation, the negative exponent were found 
among them, so the existence of the anomalous scaling was proved in the model of the 
vector advection with the mixing of the operators. It was in the case of flat flows (d = 2) 
|2*T] that the anomalous exponents of the structure functions of the arbitrary order were 
calculated in the model (jlj), (jJJ), (jHJ). The anomalous exponents were obtained in the 
one-loop approximation coincide with the exponents (j3J) of the scalar admixture. "The 
most dangerous" operators were found to be the powers of dissipation operator, just as 
in the scalar model. 

3. The bases for scalar operators of the form (d(f) n 

As demonstrated in |18j . the leading terms in (J2J) arise from the family $( n ) of scalar 
operators of the form (d(f) n . More precisely, that family consists of all possible 
contractions of n tensors d^j. Denoting the symmetric and the antisymmetric parts of 
the ditpj as 

Sy = (diifj + dj<pi)/2, = (di<pj - dj(fi)/2, (9) 

we can express for 64 the seven operators which forms the family <3>( 4 ) = {tr(A 2 ) 2 , 
ti(S 2 ) 2 , tr{A 2 )tr{S 2 ), tr{A 2 S 2 ), tr((AS) 2 ), tr(A 4 ), tr(S 4 )}. The main difficulty in the 
calculation of the critical dimensions of the operators from $( n ) is the rapid growth 
with n of the number of relevant operators. For example, the family $^ 6 - ) contains 24 
operators, $ ( - 8 ^ — 81, $( 10 ) — 278, and $( 18 ) — as many as 47246 operators. Fortunately, 
not all of them are independent for sufficiently small space dimensions d. For the most 
interesting cases d = 2 and d = 3, the elimination of redundant operators reduces 
drastically the size of the matrix of critical dimensions. 

Our goal is to construct for d = 2 and d = 3 the sets of minimal size which 
contains invariants of the matrix diifj. Then we will form bases of & n ' using products 
of the invariants from f^. 

3.1. 2-d case 

The dependencies mentioned above can be considered as consequences of the Hamilton- 
Cayley identity. For a traceless 2x2 matrix M it has the form: 

Xi{M) = M 2 — Jtr(M 2 )/2 = 0, (10) 

with identity matrix I. From ()10|) one has \2{ctA + (3S) = for any linear combination 
of the matrices A and S from @. Collecting the coefficients of three independent 
structures a 2 , (3 2 and a(3 gives three relations: 

A 2 - Jtr(yl 2 )/2 = 0, S 2 -Itr{S 2 )/2 = 0, AS + SA = 0. (11) 
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The first and the second relations are the Hamilton-Cayley identities for the matrices 
A and S. The last one can be interpreted as a commutation relation. The relations (fTTj) 
can be written in the form 

A 2 ~ 0, S 2 ~ 0, SA ~ -AS, (12) 

where ~ means the equality up to polynomials of lesser degree. 

It follows from (|T2~j) that P 3 (A, S) ~ for any polynomial of the third degree, so 
that an arbitrary polynomial of the matrices A and S can be written as P(A, S) = 
Ci(f!,2)AS + 6*2(02) A + 63(^2) S + 64(02)/, where the coefficients are polynomials of 
the scalar invariants from the minimal set 

Q 2 (A,S) = {tr(A 2 )MS 2 )} (13) 

and, as a consequence, tr(P(A, S)) = 2C4(0 2 ). 

Thus for d = 2 each operator from $( 2n ) is a polynomial P n (tr(A 2 ), tr(S' 2 )), which 
can be decomposed in the basis 

tr(A 2 ) n ~ k ti(S 2 ) k , 0<k<n. (14) 



3.2. 3-d case 

The Hamilton-Cayley identity for a traceless 3x3 matrix M has the form: 

X 3 (M) = M 3 — M tr(M 2 )/2 - / tr(M 3 ) /3 = 0. (15) 

From the Hamilton-Cayley identity for the matrices A and S from Q it follows that 
the minimal set ^(A, S) contains the invariants tr(A 2 ), tr(S' 2 ) and tr(S* 3 ). Equating to 
zero the coefficients of a 2 [3 and af3 2 in the identity Xz{ a A + (3S) = gives: 

A 2 S + ASA + SA 2 - Str{A 2 )/2 - Itr(A 2 S) = 0. (16) 

AS 2 + SAS + S 2 A - A tr(S 2 )/2 = 0. (17) 

It follows from (jlfij) that the invariant tr(A 2 S) is in the minimal set flz(A, S). Neglecting 
the polynomials of lesser degrees in (|15 p — (|17| ) results in commutation-like relations 

A 3 ~ 0, S 3 ~ 0, (18) 

S*yL4 ~ - ASA, SSA ~ - SAS. (19) 

It follows from ()18|) that each monomial containing A 3 or S" 3 vanishes. Each monomial 
containing SAA or SSA can be decomposed according to (fTHj) as a sum of other 
monomials each of which precedes the initial monomial in the lexicographic ordering. It 
then follows that the iterations of ()18j) and p9j1 finally give an expression without any 
one of the factors A 3 , S 3 , SA 2 and S 2 A. Applying (JTHJ) and (JTHJ) to a polynomial of the 
third degree gives: 

P 3 (A, S) ~ L(A 2 S, ASA, AS 2 , SAS), (20) 

where L is a linear combination of the arguments. Moreover, from relation (|2U|) it follows 
that tr{P 3 {A,S)) = L{ti{A 2 ),tr{A 2 S),tr{S 2 ),tY{S 3 )). 
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Applying (fTS|) and ()19|1 to a polynomial of the fourth degree gives 

P A {A, S) ~ L(A 2 SA, A 2 S 2 , AS AS, SASA, SAS 2 ), (21) 

but the expression allows further simplification. The reason is the existence of another 
relations besides (fT5 j) - (fT7]) . For example, none of the relations (JTHj). fjTT?^ can be applied 
to the monomial AS AS AS, but from the Hamilton-Cayley identity for the matrix AS 
it follows that (AS) 3 ~ 0. Applying (fT5j) - (fT7]) to the coefficient of a(3^ in the identity 
Xz(oiA + (3S + J AS) = gives the equation 

A 2 S 2 -(SA) 2 -A 2 tr(S 2 )/2-S 2 tr(A 2 )/2+Itr(A 2 S 2 )+Itr((AS) 2 ) = 0.(22) 

Neglecting the polynomials of lesser degrees results in 

SASA ~ A 2 S 2 (23) 

that turns flSJ into P 4 (A, 5) ~ L(A 2 SA, A 2 S 2 , AS AS, SAS 2 ). Taking the trace of p2|) 
gives the equation that relates tr(y4 2 5* 2 ) and tr((AS) 2 ): 

4 tr(A 2 5 2 ) + 2 tr((AS) 2 ) - tr(A 2 ) tr(5 2 ) = 0. (24) 

Thus we have found all possible relations for the fourth-degree monomials. From 
f!22j) and (|2*4"j) the necessity follows to add one more invariant to f^A, S) — tr(A 2 S 2 ), 
tr((AS) 2 ) or their linear combination independent on (|2*4*j) . for example, the invariant 
tr([A, S} 2 ) with the commutator [A, S] = AS — SA. 

For the polynomials of the fifth degree, applying of (fT£|) . (fTTiJ) and (f2~3|) gives 
P§(A, S) ~ L(A 2 S AS , AS AS 2 ) . To find relations for the fifth degree monomials we 
equate to zero the coefficient of a(3 2 in the identity x^( a S+[3AS) = after simplification 
that results in I ti(ASAS 2 ) + Itr(A 2 ) tr(S' 3 )/6 = 0. Consequently, there are no new 
invariants to add to the fl 3 (A, S). 

Using (|18J). (|19J) and (}2*3*j) for the polynomials of the sixth degree gives the relation 

P 6 (A,S) ~ L(A 2 SAS 2 ), specifically, it follows [A,S] 3 6A 2 SAS 2 . Analyzing the 

identity X3([A<S']) = we conclude that [A, S] 3 ~ 0. Another invariant tr(L4, S] 3 ) = 
-6ti(A 2 SAS 2 ) is added to the set fi 3 . 

Finally, P&(A, S) ~ 0, and hence any polynomial of the matrices A and S can be 
represented as a polynomial of the degree not higher than five with coefficients which are 
polynomials of f2 3 (A, S). For the trace of each polynomial we have tr(P(A, S)) = P(Q%). 
The complete minimal set of invariants of the traceless 3x3 matrix (JHJ) is 

n 3 (A, S) = {ti(A 2 ),tT(S 2 ),tT(A 2 S),tT(S 3 ),tr([A, S] 2 ),tr([A, S] 3 )}. (25) 

To construct the basis of the family $ < - n - ) we need to take into account additional relation 
among invariants in (|23|) . Namely, that the last invariant is similar to pseudoscalar and 
its square is expressed as the polynomial on the others invariants. To prove that we 
rewrite (JTHJ) in the form Jtr(M 3 )/3 = M 3 -Mtr(M 2 )/2 substituting M = [A, S]. After 
squaring and applying of the (fTT)|) . (fTTj) and (j22J) to the right hand side we obtain the 
desired relation. Thus for d = 3 any operator from $^ n ^ can be decomposed in the basis 

ti(A 2 ) ni tr(S 2 ) n2 tr(A 2 S) ns tr(S 3 )™ 4 ti([A, S} 2 ) 115 tr(L4, S} 3 ) ne , 

n = 2ni + 2n 2 + 3n 3 + 3n 4 + 4n 5 + 6n 6 , n G < 1. 
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4. Conclusion 

The leading contribution to the inertial-range behavior of the structure functions of 
the vector admixture tp has the form S n (r) ~ r n ( 1 - £ / 2 )( r/ /£)-?™ ) where the anomalous 
exponents £ n are given by eigenvalues of the matrix of critical dimensions of the family 
$("•) f scalar operators of the form (dtp) n , see [18J. The main results of the paper 
are finding of the minimal sets of the invariants of the tensor d^j for the spatial 
dimensions d = 2 (|T3|) and d — 3 (|25|) . The method proposed for invariants searching is 
based on the Hamilton-Cayley identity and can be generalized to searching of the joint 
invariants of the several tensors for arbitrary spatial dimension. Using the invariants 
from the bases of the family $( n - ) were constructed for arbitrary n in two dimensions 
(114)1 and three dimensions (|26j). 

For d — 2, the basis (fT4*)) has clear physical sense as tr(S 2 ) oc $dis - - the 
dissipation operator and tr(v4 2 ) oc |rot(/?| 2 - - the enstrophy operator. The matrix 
of the critical dimensions become triangular in the basis. Just as in the scalar 
admixture case the inertial range behaviour of the structure function is determined 
by the power of the dissipation operator. So the basis ()14)) provides the possibility to 
solve problem in general. With another representation of the basis (fH^l it was obtained 
£ n = en(n — 2)/8 + 0(e 2 ) in |2T] that coincides with the exponents (0) of the scalar case. 

For d = 3, the matrix of critical dimensions in the basis (|2fi)l has a general form. 
However, taking into account the dependencies between the operators radically reduces 
the size of the matrix. For example, the family $( 18 ) contains 47246 operators for general 
d, but only 154 of them appear independent for d = 3. This gives the possibility to 
calculate the anomalous exponents (in the linear approximation on e): £ 4 « 0.546e, 
£ 6 ~ 1.75c, & « 3.66e, £io ~ 6.27e, £ 12 ~ 9.58^, ^4 w 13.6e, 6e ~ 18.3e, 6s « 23.7^. 
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